Streamflow simulation in watersheds is necessary to runoff management and understanding watershed response. Different hydrological models were developed for rainfall-runoff modelling. The BILAN model was used to simulate daily runoff in 23 river gauge stations in a 20-year period using observed rainfall, temperature, humidity and discharge time series. The values of model parameters including Spa, Dgm, Dgw. Alf, Mec, Wic, Soc, and Gid have been calibrated using the genetic algorithm. The monthly values of streamflow data were simulated in the selected watersheds. The efficiency of the calibrated model was evaluated using statistical measures. The statistical measures were coefficient of efficiency (CE), root mean square error (RMSE) and coefficient of determination (R2). The results showed that in the calibration period, the ranges of values of R2 evaluation indicators were acceptable in Shams Abad (92%), Neneh Karan (42%), Mashiran (38%) and Barouq (51%) watersheds. Also, the results of the BILAN model validation showed that the efficiency of the model in Amoghin, Atashgah, Nir, Namin, Samian, Viladaragh and Shams Abad watersheds had a good agreement with the R2 values more than 60%. Meanwhile, the performance of the BILAN model in Gilandeh, Iril, Ahmadkandi, Barouq and Poleh Almas stations was not acceptable according to the R2 values less than 5%. In general, according to the results of model performance evaluation, it can be said that the BILAN model has a relatively acceptable performance in some of studied stations. Therefor the BILAN model can be implemented based on the available input data of the model.